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FINAL  REPORT  ON  SAI  SUPPORT  OF  THE 
ONR  UPPER  CRUST  STUDY 

The  SAI  support  of  the  ONR  Upper  Crust  Study  con- 
sisted of  two  efforts:  (1)  An  SAI  representative  (J.  S. 
Hanna)  attended  the  meetings  of  the  Crust  Study  Working 
Group  during  which  the  relevant  data  from  past  exercises 
were  reviewed  and  plans  were  made  for  a comprehensive 
geophysical/acoustic  experiment  and  (2)  Two  tutorials  on 
matters  related  to  bottom-interacting  acoustic  energy  were 
prepared  in  support  of  the  group's  planning. 

The  documentation  of  these  efforts  has  been  pre- 
pared separately  and  consists  of  the  following: 

(1)  Draft  Final  Report  on  the  Discussions  of  the 
Ocean  Crust  and  Lithosphere  Working  Group 

(2)  Draft  Proposal  for  the  Rivera  Ocean  Seismic 
Experiment  (ROSE) 

(3)  "Reflection  Loss  for  Long-Range  Transmission 
Loss  Estimates:  What  Matters  and  How  it  Might 
be  Measured,"  SAI  Report  No.  76-688-WA 

(4)  "Some  Complications  in  the  Traditional  Mea- 
surements of  Bottom  Loss,"  SAI  Report  No. 
76-664-WA 

Items  1 and  2 above  have  been  assembled  and  distributed  by 
J.  Ewing  and  are  based  upon  the  efforts  of  the  Working  Group. 
Items  3 and  4 were  prepared  by  J.  S.  Hanna  of  SAI  in  partial 
fulfillment  of  Contract  N00014-76-C-1049  and  are  enclosed 
with  this  report. 
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INTRODUCTION 


In  the  past  few  years,  two  substantial  measurement 
programs^  ’“  have  been  conducted  for  the  purpose  of  measuring 
bottom  loss  in  the  low-frequency  range  below  1 kHz.  Both 
programs  have  processed  signals  from  explosive  sources  in 
various  one-third  octave  bands  and  have  used  similar  source- 
receiver  geometries  with  both  source  and  receiver  within 
300  to  800  feet  of  the  ocean  surface.  For  frequencies  below 
100  Hz.  it  is  shown  here  that  significant  artifacts  may  be 
induced  in  the  inferred  bottom  loss  because  of  the  combined 
choice  of  signal  processing  and  source-receiver  geometry. 
Furthermore,  over  the  entire  frequency  range  of  interest, 
it  is  often  likely  that  the  inferred  losses  at  low  grazing- 
angle  are  contaminated  by  unsuspected  intrusions  of  energy 
from  the  first  convergence  zone. 
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MEASUREMENT  BACKGROUND 


/ 
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It  has  become  widely  appreciated  ' ' that  at  some 
low  frequency  (probably  below  1 kHz),  it  becomes  inappropriate 
to  consider  that  all  energy  incident  upon  the  ocean  bottom 
is  simply  reflected  at  this  interface;  rather,  one  must 
account  for  the  energy,  transmitted  through  this  boundary, 
which  is  returned  to  the  water  column  by  reflection  or  re- 
fraction in  the  ocean  bottom.  In  fact,  at  frequencies  of  a 
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few  hundred  Hertz  and  below  for  the  unconsolidated  sedi- 
mentary bottoms,  limited  data  suggest  that  relatively  little 
energy  is  reflected  from  the  small  impedance  discontinuity 
at  the  bottom  interface  and  the  most  important  energy  paths 
are  reflected  and  refracted  in  the  bottom.  It  is  possible 
that  the  attenuation  of  sound  in  the  ocean  bottom  makes 
these  sub-bottom  paths  significant  onlv  for  these  low  fre- 
ouencies. 


To  preserve  some  generality  in  the  discussion  of 
principles  which  follows,  then,  a sound-speed  profile  which 
Includes  the  sediment  layer  below  the  ocean  bottom  (Figure  1) 
is  assumed.  The  sound-speed  in  the  water  column  is  typical 
of  the  north  Pacific,  while  that  in  the  sediment  is  character- 
ized bv  a constant  gradient  of  1 sec-1  which  is  representative 
of  measurements  for  sediments.  The  perturbing  influences 
of  such  quantities  as  density  discontinuities  and  attenuation 
in  the  sediment  will  be  included  in  a qualitative  way  where 
appropriate . 

The  measurements  of  References  1 and  2 are  actually 
of  transmisison  loss;  a nrona^ar ton-loss  model  is  then  assumed 
for  the  bot tom- interact ing  paths  involved  and  the  difference 
between  measured  and  computed  loss  is  imputed  to  reflection 
loss  at  the  bottom.  The  measurements  attempt  to  isolate 
the  four  principal  oaths  shown  schematically  in  Figure  2. 

All  the  examples  discussed  here  will  be  for  a source  depth 
of  800  feet  and  receiver  depth  of  300  feet.  On  the  scale 
employed  in  Figure  2.  it  would  not  be  possible  to  distinguish 
graphically  whether  the  paths  reflect  from  the  bottom  or 
are  refracted  by  the  assumed  strong  gradient  in  the  sediment. 
In  the  following  discussion,  however,  the  distinction  will 
be  carefully  maintained. 
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3.0  COMPLICATIONS  CAUSED  BY  THE  GEOMETRY  AND 

SIGNAL  PROCESSING 

To  gain  a first-order  appreciation  for  the  con- 
sequences of  geometry  and  signal  processing,  consider  the 
four  paths  of  Figure  2 and  assume  that  they  are  refracted 
in  the  sediment  (rather  than  reflected  from  the  bottom). 

For  the  sources  typically  used  (1.8  lb  Mk  61  SUS  charges 
detonated  at  800  feet),  the  duration  of  the  source  wave- 
form does  not  permit  resolving  these  paths,  in  general, 
and  so  both  measurement  programs  have  assumed  that  these 
paths  contribute  (1)  equally  and  (2)  incoherently  to  the 
measured  loss.  The  first  assumption  is  not  unreasonable; 
however,  the  discussion  below  demonstrates  that  for  the  geo- 
metries employed  and  for  frequencies  of  100  Hz  or  less, 
the  second  assumption  is  inappropriate. 

The  key  to  understanding  the  inadequacy  of  the 
second  assumption  is  the  relative  travel  times  for  the  four 
paths.  These  are  plotted  in  Figure  3 for  the  range  interval 
from  14  to  29.5  nautical  miles  which  covers  the  bottom  grazing 
angles  from  20  down  to  2 degrees  for  the  assumed  sound-speed 
of  Figure  1 and  geometry  of  Figure  2.  Since  the  measurement 
consists,  at  any  range,  of  an  average  over  frequency  of  the 
impulse  response  associated  with  this  four-path  structure, 
it  is  instructive  to  consider  several  examples  of  the 
impulse  response. 

Assuming  the  paths  to  be  of  equal  amplitude  and  the 
relative  phases  to  be  determined  only  by  the  relative  travel- 
times  of  Figure  3 and  the  phase  reversals  from  surface  re- 
flections, the  Fourier  transform  of  the  impulse  response  is 
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where  u is  the  circular  frequency  and  T0,  Tg,  and  are 
the  relative  travel-times  shown  in  Figure  3.  Notice  in 
Figure  3 that 

T 2 “ T4“T3  and  T3=T4-T2 

Thus,  there  are  really  only  two  basic  travel-time  differences, 
the  shorter  corresponding  to  the  up-  and  down-going  paths 
at  300  feet  and  the  longer  to  the  up-  and  down-going  paths 
at  800  feet.  Plots  of  I l|“  versus  frequency  will  thus  show 
periodicities  of  1/T0  and  1/T^.  Figures  4 and  5 show  such 
plots  of  20  log  |l|  for  the  extreme  ranges  of  interest,  14 
and  29.5  nautical  miles. 

To  appreciate  the  inadequacy  of  assuming  that  the 
paths  add  incoherently,  imagine  averaging  | 1 1 2 in  either 
Figure  4 or  5 over  one-third  octave  intervals.  Such  in- 
tervals vary  in  width  from  about  9 Hz  at  a center  frequency 
of  35  Hz  to  23  Hz  at  a center  frequency  of  100  Hz.  The 
assumption  of  incoherent  addition  is  equivalent  to  assuming 
that  the  one-third  octave  average  of  |l|  yields  approxi- 
mately 1/4.  Figures  4 and  5 make  it  clear  that  this  will  not 
be  the  case,  in  general,  for  the  geometry  discussed  here 
and  the  frequency  interval  shown. 


Figure  6 illustrates  the  effect  on  inferred  re- 
flectivity of  the  partially  coherent  summation  of  paths 
discussed  above.  For  the  four  paths  under  consideration, 
the  bottom  panel  of  Figure  6 shows  computed  total  transmis- 
sion loss  assuming  (1)  incoherent  (RMS)  summation,  (2)  a 
proper  one-third  octave  average  at  35  Hz  and  (3)  a one- 
third  octave  average  at  100  Hz.  If  we  assume  that  the  model 
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of  four  paths  embodies  all  that  matters  in  a hypothetical 
measurement  and  that  the  one-third  octave  averages  corre- 
spond to  data,  the  top  panel  of  Figure  <5  shows  the  inferred 
reflectivity  at  35  and  100  He;  using  the  KMS  sum  for  the 
the  assumed,  or  modeled,  transmission  loss.  There  is.  of 
course,  no  bottom  loss  in  the  "data';  however,  this  figure 
displays  the  spurious  structure  induced  in  the  inferred 
reflectivity  by  the  failure  to  model  accurately  the  summa- 
t toil  of  paths  appropriate  to  the  geometry  and  signal  pro- 
cessing employed.  Ad  the  frequency  decreases,  the  spurious 
structure  varies  slower  with  angle  and  is  more  exaggerated. 


That  the  effect  discussed  above  exists  in  the  mea- 
surements is  suggested  by  the  data  in  Figure  7.  These  data 
are  from  a Naval  Oceanographic  Office  exercise;  the  fre- 
quency was  50  Hz , the  source  depth  was  $00  feet,  and  the 
receiver  depth  was  300  feet.  For  comparison,  a computation 
for  the  parameters  of  the  measurement  (.like  those  in  Figure  0) 
is  also  shown.  Clearly,  the  data  show  essentially  the  same 
structure  as  the  computation.  The  last  two  data  points 
(beyond  27  nml  are  felt  to  be  evidence  of  another  complica- 
tion which  is  discussed  in  the  next  section. 

The  discussion  thus  far  has  assumed  that  all  energy 
incident  upon  the  water-bottom  interface  was  transmitted 
through  the  boundary  and  refracted  by  the  gradient  back  into 
the  water.  In  principle,  however,  partial  reflections  will 
occur  at  the  boundary.  This  additional  mechanism  would 
result  in  eight,  rather  than  four,  paths  with  additional  time- 
delays.  each  causing  additional  interference  patterns.  It  is 
expected,  then,  that  measurements  may  reveal  more  structure 
than  the  four-path  computation  This  expectation  seems  born 
out  by  the  data  of  Figure  8.  These  data  are  from  the  same 
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exercise  as  those  of  Figure  7 and  correspond  to  a frequency 

I of  100  Hz.  These  data  seem  to  contain  not  only  the  structure 

of  the  four-path  computation,  but  more  besides.  Again,  the 
, last  two  data  points  will  be  interpreted  as  evidence  for 

an  effect  discussed  in  the  next  section. 

T 
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4.0  COMPLICATIONS  CAUSED  BY  WATERBORNE  ENERGY 


I 

i 

1 

I 

1 


Up  to  this  point,  it  has  been  assumed  that  the  four 
bottom- interact ing  paths  can  be  resolved  from  all  other  paths 
in  the  problem.  This  is  not  always  the  case  for  ranges  cor- 
responding to  low  grazing  angles.  To  demonstrate  this  fact, 
consider  first  the  ray  plot  of  Figure  9 where  are  shown  the 
rays  from  a source  at  800  feet  which  arrive  in  the  range- 
depth  window  from  25  to  35  nm  and  0 to  300  meters.  Those 
paths  which  reflect  from  the  surface  are  dist inguished 
according  to  whether  they  interact  with  the  bottom  or  corre- 
spond to  the  set  of  waterborne  rays.  (The  notation  RSR 
stands  for  refracted/surface-reflected  and  denotes  paths 
which  reflect  from  the  surface  and  turn  at  depths  above  the 
ocean  bottom. ) Consider  the  rays  which  intersect  the  re- 
ceiver depth  at  300  feet:  the  last  bottom-interacting  path 
arrives  at  a range  of  31  nm,  yet  even  the  ray-trace  shows 
non  bottom-interacting  paths  arriving  in  the  overlapping 
range  from  29.5  to  31  nm.  In  reality,  however,  the  water- 
borne paths  make  their  influence  felt  before  29.5  nm  in  the 
form  of  the  shadow  zone  field  of  the  caustic.  The  relative 
travel  time  between  the  refracted  field  and  the  bottom  inter- 
acting field  is  sufficiently  small  so  as  not  to  be  resolved 
by  1/3  octave  processing  at  low  frequencies.  Thus,  attempts 
to  measure  the  bottom-interacting  field  at  these  ranges  may 
be  thwarted  by  the  additional  influence  of  the  refracting 
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To  estimate  this  effect,  a wave  model  which  properly 
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treats  caustic  fields  (the  parabolic  equation  model  ' ) was 
run,  including  both  the  bottom-interacting  and  refracted 
fields,  and  averaged  in  frequency  over  1/3  octave  at  35  Hz, 
The  results  of  these  calculations  are  shown  in  Figure  10. 

To  estimate  the  reflectivity  shown  at  the  top  of  the  figure, 
the  parabolic  equation  computation  is  assumed  as  data  and 
the  RMS  sum,  again,  as  estimated  loss;  the  reflectivity 
follows  the  earlier  35  Hz  curve  down  to  about  5 deg. , below 
which  it  goes  even  further  into  negative  values  than  before. 
The  significant  point  here  is  that  even  if  the  proper  summa- 
tion of  the  bottom-interacting  paths  were  used  as  the  esti- 
mated transmission  loss,  negative  reflectivity  losses  would 
be  obtained  at  low  grazing  angles  because  of  the  refracted 
contribution  to  the  field.  The  extent  of  the  refracted  field 
will  depend  upon  frequency,  geometry  and  depth  excess. 

The  50  Hz  data  in  Figure  7 are  felt  to  give  evidence 
of  the  contamination  of  the  waterborne  energy  discussed 
above.  To  illustrate,  the  data  and  calculation  of  Figure  7 
are  reproduced  in  Figure  11  with  the  addition  of  a parabolic 
equation  estimate  for  the  contribution  of  the  waterborne 
energy  alone.  Note  that  as  range  increases  to  26  nm,  the 
measured  loss  increases.  Beyond  this  range  the  loss  de- 
creases rapidly  and  follows  almost  precisely  the  estimated 
loss  for  only  waterborne  energy.  A possible  interpretation 
of  the  data  now  appears  to  be  that  (1)  the  bottom  interacting 
contribution  is  falling  off  at  a rate  suggested  by  the  data 
from  20  to  26  nm  and  (2)  the  received  field  beyond  26  nm  is 
dominated  by  the  refracted  contribution. 
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The  100  Hz  data  of  Figure  8 are  likewise  felt  to 
display  the  influence  of  this  waterborne  energy.  Those  data 
are  reproduced  in  Figure  12  along  with  estimates  for  t h» 
waterborne  energy  alone  and  the  appropriate  one-third  octave 
coherent  combination  of  both  waterborne  and  bottom  interacting 
energy.  Again,  beyond  27  nm  the  data  seem  quite  consistent 
with  the  additional  contribution  of  the  waterborne  energy. 


5.0  CONCLUSIONS 

Examination  of  measured  transmission  loss  data  which 
are  used  to  infer  bottom  reflectivity  shows  strong  evidence 
of  spurious  structure  in  the  inferred  reflectivity.  The 
origins  of  this  structure  are  found  in  (,11  the  geometry  of 
the  measurements  and  (2)  the  signal  processing  of  the  data. 

A basic  shortcoming  in  the  measurements  is  the  inability  to 
distinguish  only  that  energy  which  interacts  with  the  bottom 
and  to  resolve  the  multipath  structure  caused  by  the  shallow 
source/ receiver  geometry.  Recent  work  by  Di oils'  ’ eliminates 
the  geometry- induced  multipath  problem;  however,  any  source 
receiver  geometry  may  have  the  complication  of  waterborne 
energy  contaminating  the  measurements  for  low  grazing  angles 
of  the  bot tom- Interact ing  energy. 
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Figure  3:  Relative  Travel  Times  for  Profile  in  Figure  1 


Figure  4:  Impulse  Response  Spectrum  at  14nm 
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Figure  6:  Effect  of  Geometry  and  Processing 
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Figure  9:  Ray-Trace  Range/Depth  Hindow  of  the  First  Convergence  Zone 


E ENERGY 


REFLECTION  LOSS  FOR  LONG-RANGE 
TRANSMISSION  LOSS  ESTIMATES: 
WHAT  MATTERS  AND  HOW  IT 
MIGHT  BE  MEASURED 


ATLANTA  • ANN  ARBOR  • BOSTON  • CHICAGO  • CLEVELAND  • DENVER  • HUNTSVILLE  • LA  JOLLA 
LITTLE  ROCK  • LOS  ANGELES  • SAN  FRANCISCO  • SANTA  BARBARA  • TUSCON  • WASHINGTON 


REFLECTION  LOSS  FOR  LONG-RANGE  TRANSMISSION 
LOSS  ESTIMATES:  WHAT  MATTERS  AND  HOW  IT 
MIGHT  BE  MEASURED 


SA I- 76-688- WA 


November  1976 


Prepared  for: 

Office  of  Naval  Research 
Arlington,  Virginia  22217 


Prepared  by: 
John  S.  Hanna 


Prepared  under:  Contract  N00014-76-C-1049 


SCIENCE  APPLICATIONS,  INC. 

8400  Westpark  Drive,  McLean,  Virginia  22101 
Telephone  703/821-4300 


1.0 


INTRODUCTION 


The  consideration  of  the  acoustic  influence  of 

the  ocean  bottom  ranges  all  the  way  from  a characterization 
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of  this  interface  by  a reflection  coefficient  ’ (which  may 

be  a catchall  for  several  loss  mechanisms),  to  detailed 

attempts  at  including  the  basic  geophysics  of  the  bottom 
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in  one's  computations  ’ (e.g.,  sediment  sound  speed, 

density,  absorption,  etc.).  Many  models  which  attempt 
to  make  estimates  of  average,  rather  than  detailed,  acoustic 
propagation  use  reflection  coefficients  which  also  represent 
certain  average  properties  of  the  ocean  bottom.  In  the 
context  of  these  kinds  of  calculations  it  is  instructive 
to  understand  what  features  of  the  reflection  coefficient 
are  important  and,  from  there,  to  consider  how  they  might 
be  inferred  from  acoustic  data.  The  following  discussion 
is  devoted  to  such  an  understanding.  It  considers  only 
the  case  of  long-range  propagation. 


2.0  THE  IMPORTANCE  OF  BOTTOM  INTERACTIONS 


Existing  measurements  of  bottom  reflection  coef- 
ficients below  several  hundred  Hertz  show  quite  a wide 
range  of  loss-versus-angle  functions.  Some  measurements 
have  shown  the  loss  to  be  6 to  8 dB  at  all  angles  while 
others,  to  the  accuracy  of  the  measurement,  show  no  loss 

from  0 degrees  grazing  up  to  some  critical  angle,  with 

1 2 

a rapidly  increasing  loss  beyond.  ’ 

Just  to  demonstrate  the  implications  of  such 
variations  in  reflection  coefficient,  two  propagation  loss 

5 

estimates  were  made  using  FACT,  for  a typical  Pacific 


velocity  profile  in  deep  (18,000  feet)  water  corresponding 
to  the  reflection  coefficient  extremes  mentioned  above. 
These  computations  are  summarized  in  Figures  1 and  2;  the 
inset  in  the  upper  right  corner  of  each  figure  gives  the 
assumed  coefficient. 

In  Figure  1,  for  the  high-loss  coefficient,  the 
contribution  to  the  propagation  from  energy  which  interacts 
with  the  bottom  is  seen  to  fall  off  very  rapidly  with 
range,  leaving  only  the  contribution  from  the  waterborne 
energy  which  appears  in  the  form  of  the  usual  convergence 
zones.  By  contrast,  the  results  in  Figure  2,  for  the  low- 
loss  case,  show  that  the  bottom-interacting  energy  domin- 
ates the  propagation  to  long  range.  Both  these  extremes 
have  been  verified  in  actual  measurements  of  propagation 
loss. 


Given  such  extreme  possibilities  it  becomes  inter- 
esting to  ask  what  features  of  a reflection  coefficient  are 
most  important:  the  critical  angle,  the  functional  shape, 
or  what?  Clearly,  one  approach  would  be  to  build  a compen- 
dium of  cases  by  performing  a parametric  study.  An  alter- 
native, pursued  below,  instead  makes  some  reasonable 
simplifying  assumptions  which  permit  derivation  of  an 
expression  for  the  bottom-interacting  energy  that  is 
amenable  to  intuitively  satisfying  interpretation;  it 
also  suggests  some  interesting  interpretations  of  measured 
propagation  loss  for  the  purpose  of  estimating  the  effec- 
tive reflection  loss. 
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3.0 


THE  MODEL 


From  the  point  of  view  of  geometrical  optics, 
those  paths  which  interact  with  the  bottom  are  insensitive 
to  (1)  refraction  by  the  water  column  which  seldon  affects 
in  a significant  way  the  spreading  loss  along  the  path 
that  is  otherwise  determined  simply  by  the  path  length 
and  (2)  the  source  and  receiver  depths,  since  for  all 
choices  (at  least  in  water  which  is  not  bottom-limited), 
the  same  bottom-interacting  paths  will  connect  the  two 
depths.  Based  on  these  two  points,  the  simple  model  of 
Figure  3 is  used.  The  water  column  is  assumed  to  be 
isovelocity  so  that  the  ray  paths  are  straight  lines; 
these  straight  lines  are  reasonable  approximations  for 
bottom  interacting  paths  even  in  the  presence  of  refraction. 
The  existence  of  refraction  for  the  real  ocean  environment 
has  essentially  only  one  important  consequence  for  the 
model  of  Figure  3:  there  is  a maximum  permissible  ray 
period  for  bottom  interacting  paths  which  means  we  should 
only  consider  angles,  0,  down  to  some  minimum.  ^min . in 
the  isovelocity  model.  For  a typical  deep-ocean  case  this 
angle  is  of  order  ten  degrees.  This  restriction  is  just 
another  way  of  saying  that  one  important  quantity  is  the 
angular  aperture  at  the  source  corresponding  to  bottom 
interacting  paths. 

Before  postulating  a functional  form  for  the 
bottom  reflectivity  it  is  instructive  to  pursue  the  above 
point  a bit  further.  In  our  isovelocity  model  the  ray 
angle  (which  is  constant  along  a path)  should  be  viewed 
as  the  angle  at  the  source  in  the  real  world.  For  example, 
9min  in  our  model  w°uld  correspond  to  the  real-world  ray 
which  just  grazes  the  bottom.  When  considering  the 


correspondence  between  the  angle  6 in  the  model  and  the 
angle  at  the  bottom,  t?B,  for  the  real  world,  the  following 
transformation  is  appropriate: 


„ _ COS  0 

cos  9n  - , ■ -a 

B cos  emin 


This  transformation  is  just  an  application  of  Snell's 
Law  assuming  that  the  angle  ®min  corresponds  to  90  ■ 0. 

The  ratio,  a,  of  the  reflected  intensity  to  the 
incident  intensity  is  modeled  as 


a(0)  - exp 


' Yl(SJn  "*)' 


This  expression  is  applied  up  to  some  critical  angle, 

0 , , beyond  which  a(9)  is  assumed  to  decrease  rapidly  to 
some  negligibly  small  quantity:  y^  and  Yj  are  constants, 
In  this  expression  the  quantity 


exp  ( -yq) 


is  the  ratio,  a,  at  grazing  (.i.e.,  0 
quantity 

exp  ;Y°  - Yi  fe " 0 


is  the  ratio  at  the  critical  angle. 


0m<„)  and  the 
min 
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The  functional  form  in  Equation  2 was  chosen 
primarily  to  facilitate  the  closed  form  computation  of 
intensity  as  a function  of  range.  In  many  parameteriza- 
tions  the  reflection  loss,  in  dB,  is  assumed  to  vary 
linearly  with  angle;  such  an  assumption  would  require 
the  following  form  for  the  intensity  ratio; 

e*p  £-Y0  - Yi  <«  - 8min)J  (3) 

However,  this  form  does  not  permit  a closed  form  analysis 
of  the  type  presented  here;  furthermore,  the  transformation 
of  Equation  1 to  take  Equation  3 from  our  model  angle  0 to 
the  real-world  angle  0Q  would  result  in  an  effective  reflec- 

D 

tion  loss  (in  dB)  which  is  not  linear  in  0g.  So,  the 
interesting  question  is  whether  Equation  2,  when  trans- 
formed from  0 to  0q,  gives  a loss-versus-angle  of  some 
acceptable  form.  As  an  example,  Figure  4 shows  the  effec- 
tive loss  as  a function  of  0g  when  Equation  2 is  required 
to  produce  no  loss  at  grazing  (0  = ^min)  and  2 dB  loss  at 
0 ■ 30  degrees.  Although  the  result  is  not  a straight 
line  between  the  specified  points,  it  is  deemed  acceptable 
for  the  general  questions  of  interest  here. 

In  the  Appendix  it  is  shown  that  the  model  of 
Figure  3,  along  with  the  assumed  intensity  ratio  of  Equa- 
tion 2,  yields  the  following  expression  for  intensity  as 
a function  of  range  (for  R >>  D). 
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where 


r 


I 

1 

I 

i 

i 

i 

i 

i 

] 

1 


L 


1 


T 

I 


1 

1 

1 


Ad 


e 


c 


^min ' 


D - water  depth, 
R ■ range 


Several  interesting  points  can  be  made  from  Equation  4. 
First,  if  there  is  no  loss  over  the  angular  interval  from 
dmin  to  8,  (i.e.,  y - 0),  then  the  appropriate  limit  of 
Equation  4 is 


2A6 

RD 


This  expression  says  that  the  intensity  of  the  bottom- 
interacting energy  falls  off  oylindrically  with  range  and 
is  proportional  to  the  angular  aperture,  Ad.  at  the  source 
which  is  a measure  of  energy  radiated  into  these  paths. 
Second,  if  there  is  no  loss  at  grazing  (y^  * 0),  then 
Equation  4 approaches  the  following  limit  for  large  R: 


4d 


I - 


mxn 


R Y 


which  says  the  intensity  falls  off  spherically;  this  is 
consistent  with  losing  all  but  the  most  grazing  ray  at 
long  range.  Finally,  for  constant  nonzero  loss  at  all 
angles  Equation  4 approaches  the  limit 
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4 exp  (-  ^jr^) 

I - ^ 

R“y 

which  is  also  consistent  with  losing  all  but  the  shallowest 
ray  whose  geometrical  spreading  loss  (1/R:)  is  modified  by 
the  exponential  loss.  This  exponential  factor  arises  from 
the  loss  of  a fixed  percentage  of  energy,  exp  (-Yq).  for 
every  interaction  with  the  bottom. 

Figure  5 gives  some  numerical  results  for  a 
typical  case  using  Equation  4.  (In  both  Figures  5 and  6. 
the  reflection  loss  shown  is  modeled  using  Equation  2 by 
matching  the  loss  at  grazing  and  at  critical  angles.  ) The 
critical  angle  is  assumed  to  be  20  degrees  and  the  water 
depth.  2.5  nm.  The  figure  contains  an  example  for  each 
of  the  special  cases  cited  in  Equations  5 through  7 above. 
Note  the  pronounced  sensitivity  of  these  curves  to  the 
presence  or  absense  of  any  loss  over  the  interval  from 
grazing  to  critical  angles  for  ranges  of  a few  hundred 
miles  or  more.  We  see  in  this  figure  the  first  result  of 
importance  to  reflection  coefficient  measurements:  if 
this  quantity  must  be  known  to  a fraction  of  a decibel, 
then  measurements  of  it  for  just  one,  or  even  a few. 
bottom  interactions  cannot  possibly  yield  these  accuracies. 

To  illustrate  the  relative  sensitivity  of  bottom- 
reflected  energy  to  the  other  parameter  (viz.,  critical 
angle),  the  results  of  Figure  6 are  presented.  Here  are 
plotted  the  cases  for  critical  angles  of  20  and  40  degrees 
and  constant  (with  angle)  losses  of  0 and  0.5  dB.  Note 
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that  at  a range  of  400  nm  it  is  at  least  as  important  to 
know  the  loss  per  interaction  with  the  bottom  as  it  is  to 
know  where,  between  20  and  40  degrees,  the  critical  angle 
lies. 


4.0  SOME  THOUGHTS  ON  REFLECTIVITY  MEASUREMENTS 

The  previous  section  demonstrated  that  the 
accuracies  required  for  reflection  loss  measurements  used 
in  long-range  transmission  loss  calculations  exceed  the 
achievable  accuracies  with  measurements  over  a few  inter- 
actions. This  suggests  that  only  analysis  of  long-range 
transmission-loss  data  themselves  will  yield  effective 
reflection  losses  of  required  accuracy.  In  principle  one 
could  compare  measured  transmission  loss  for  bottom  inter- 
acting energy  with  curves  such  as  those  in  Figures  5 and 
6 to  pick  effective  critical  angles  and  losses.  A prac- 
tical difficulty  may  stand  in  the  way  of  this  procedure, 
however.  It  often  happens  that  one  ill-specified  quantity 
in  measurements  of  transmission  loss  using  explosives 
(which  most  long  range  measurements  use}  is  the  source 
level.  Changing  this  level  simply  shifts  the  measured 
curve  up  and  down.  For  curves  like  those  shown  in  Figures 
5 and  6 this  could  confuse  the  determination  of  the  effec- 
tive critical  angle  since  the  critical  angle  affects  the 
calculations  in  a like  manner,  at  least  for  the  first  few 
hundred  miles.  The  shapes  of  the  curves  with  range  can. 
however,  be  compared  without  regard  for  this  problem  and 
might  yield  a useful  estimate  of  the  loss  per  interaction 
perhaps  even  as  a function  of  angle. 
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If  the  shapes  have  been  matched  to  get  this  loss, 
then  it  may  be  possible  to  derive  an  estimate  of  the  crit- 
ical angle  from  reflection  loss  measurements  made  over  one 
(or  a few)  interactions.  This  process  should  involve  only 
analysis  of  the  shape  of  the  reflection  loss  curve  without 
regard  for  absolute  values.  This  assumes  that  the  reflection 
loss  exhibits  a critical  angle,  which  is  loosely  character- 
ized as  an  angle  at  which  the  slope  of  the  reflection  loss 
increases  markedly.  Depending  upon  the  geometry  of  these 
reflection  loss  measurements  (i.e.,  whether  thay  are  self- 
calibrating or  not)  these  absolute  levels  may  also  depend 
on  accurate  knowledge  of  the  source  levels.  Thus,  it 
appears  possible  to  combine  the  results  of  two  measurements, 
each  of  which  may  be  sensitive  in  its  own  way  to  source 
level  errors,  to  get  usefully  accurate  measures  of  both  the 
critical  angle  and  the  loss  as  a function  of  angle. 
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APPENDIX 


In  the  model  of  Figure  3,  it  is  implicit  that  at 
long  range  (R  >>  D)  the  energy  associated  with  bottom- 
reflected  paths  is  uniformly  distributed  in  depth.  Without 
loss  of  generaltiy,  then,  we  will  assume  that  the  source 
and  receiver  are  at  the  same  depth  just  to  simplify  the 
bookkeeping.  If  S(9)  is  the  path  length  between  successive 
surface  and  bottom  reflections,  then  the  total  path  length 


S(0)  = 2nBS(9) 


where  n-  is  the  number  of  bottom  reflections.  Furthermore, 


we  have 


R = S(0)  cos  0 


The  intensity  along  the  ray  is 


1 cos  0 

2 — 


(A-l) 


We  want  to  sum  over  all  paths  at  some  range,  R, 
for  angles  greater  than  6min-  First  we  need  to  estimate 
the  angular  separation  between  successive  orders.  The 
geometry  of  our  problem  yields 


tan  0 * -g  • 2ng 


(A-2) 


A-l 


Taking  differentials 


A tan  0 


A0 


cos 


+ 1) 


or 


A0  * 2 | cos2 
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for  the  angular  separation  of  successive  orders.  Thus, 
the  number  of  arrivals  at  range  R per  unit  angle  is 


1_ 

A9 


2D  cos  9 


(A-3) 


If  we  now  approximate  the  following  sum  over 

arrivals 


N 


2 

cos  0 ^ 

35 

R 


by  an  integral  using  the  weighting  of  Equation  A-3,  we 
have 
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2D  cos  9 


d0 


(A — 1 ) 


A-2 


Before  simply  performing  the  integral  of  Equation 
A-4  it  is  necessary  to  modify  it  to  include  reflection  loss 
For  the  moment,  let  a be  the  ratio  of  reflected  to  incident 
intensity.  For  order  n0  the  total  intensity  diminution 
from  bottom  reflections  is  a11®.  From  Equation  A-2  we  have 

nB  " h tan  0 

so  that  as  a function  of  angle 
anB  „ tan  6/2D 


To  do  the  following  integrals  in  closed  form  it  will  be 
further  useful  to  assume  that  our  angles  of  interest  permit 
replacing  tan  0 by  0: 

nn  R0/2D 

a B “ a ' 


Inserting  this  result  into  Equation  A-4  gives 
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/max 

01 


R0./ 2D 


d0 


min 


As  stated  in  the  text,  we  will  assume  a to  be 
of  the  form 


exp 


-Y° ' y>  fe  ’ *) 


A-3 


Substituting  this  assumption  into  Equation  A-5  we  have 
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, r max  . . 

1 a 2Sd  exp  ( 25)  J exp  - (y0  + ei  M de 

9min 

. exP  ( Sr)  -2D 


!RD  R<Y0  + V6mi„ 


[-<Y0  + Y1 


/9min> 


9 . 

1 min 


Straightforward  substitution  of  the  limits  and  rearrange- 
ment of  resulting  factors  yields 


( Is! 

\"  21 


1 - exp 


yRA9 

2D 


where 


A9  “ 9 - 9 

max  min 


Y * Y0  + V9min 


Since  only  one  of  the  four  possible  paths  of  each  order 
connecting  source  and  receiver  has  been  included  in  the 
computations  thus  far,  the  total  intensity  will  be 
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